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Abstract 

A theory of spin-polarized electron transport in ferromagnet-semiconductor heterostruc- 
tures, based on a unified semiclassical description of ballistic and diffusive transport in 
semiconductors, is outlined. The aim is to provide a framework for studying the inter- 
play of spin relaxation and transport mechanism in spintronic devices. Transport inside 
the (nondegenerate) semiconductor is described in terms of a thermoballistic current, in 
which electrons move ballistically in the electric field arising from internal and external 
electrostatic potentials, and are thermalized at randomly distributed equilibration points. 
Spin relaxation is allowed to take place during the ballistic motion. For arbitrary poten- 
tial profile and arbitrary values of the momentum and spin relaxation lengths, an integral 
equation for a spin transport function determining the spin polarization in the semicon- 
ductor is derived. For field-driven transport in a homogeneous semiconductor, the integral 
equation can be converted into a second-order differential equation that generalizes the 
spin drift-diffusion equation. The spin polarization in ferromagnet-semiconductor struc- 
tures is obtained by matching the spin-resolved chemical potentials at the interfaces, with 
allowance for spin-selective interface resistances. Illustrative examples are considered. 
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1. Introduction 

In spintronics research, particular emphasis is currently placed on the study of spin- 
polarized electron transport in heterostructures formed of a nonmagnetic semiconductor 
and two (metallic or semiconducting) ferromagnetic contacts [1, 2, 3, 4]. For the actual 
design of spintronic devices, a detailed theoretical understanding of this kind of transport 
problem is indispensable. Several pertinent studies have been performed so far, which 
mostly rely on the drift-diffusion model. A number of important results have emerged, 
(i) For an interface between a metallic ferromagnet and a semiconductor without spin- 
selective interface resistance, the injected current spin polarization is predicted to be very 
low owing to the large conductivity mismatch [5]. (ii) Spin-selective interface resistances 
arising from tunnel barriers or Schottky barriers can greatly enhance the injection efficiency 
[6, 7, 8, 9, 10, 11]. (iii) A similar effect is to be expected when a sufficiently high electric field 
is applied across the semiconductor [12]. (iv) Under conditions where, in the semiconductor, 
ballistic transport prevails over drift-diffusion, the injection efficiency is controlled by the 
Sharvin interface resistance [13] unless spin-selective interface resistances are introduced. 

While the theory of spin-polarized electron transport in ferromagnet-semiconductor 
heterostructures has reached a level of considerable sophistication, it appears that certain 
aspects of the semiconductor part of the transport problem require a more systematic, 
unified treatment, such as the interplay of spin relaxation and transport mechanism all the 
way from the diffusive to the baUistic regime, and the effect of the detailed shape of the 
electrostatic potential profile. In this paper, we outline the principal ideas of a theory that 
meets these requirements. For illustrative purposes, we present calculated results for the 
position dependence of the zero-bias current spin polarization along a heterostructure as 
well as for the injected polarization as a function of bias. A detailed account of the formal 
development as well as specific applications of our theory will be published elsewhere [14] . 

2. Thermoballistic current 

Our treatment of spin-polarized electron transport relies on our previously formulated 
unified semiclassical description of (spinless) electron transport in parallel-plane semicon- 
ductor structures [15]. The description is based on the concept of a "thermoballistic elec- 
tron current" which combines elements of the baUistic and diffusive transport mechanisms. 
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Here, we briefly summarize this concept. 

Assuming a one-dimensional geometry, we consider a (nondegenerate) semiconducting 
sample enclosed between two plane-parallel ferromagnetic contacts at x — xi and x — X2, 
respectively, so that S = X2 — xi is the sample length (sec Fig. 1). The electron current 
density J{x',x") across the "ballistic interval" [a;', a;"] between two equilibration points x' 
and x" is given by 



J{x',x")=v,N,e-P'^^^^'^^"^ 



{xi < x' < x" < X2), which is the difference of the ballistic current injected into the interval 
at its left end at x' and the analogous current injected at its right end at x". The function 
(x{x) is the chemical potential at the equilibration point x. Furthermore, = {2nm* fj)^^^'^ 
is the emission velocity, iVg = 2(27rm*//3/i^)^/^ is the effective density of states at the 
conduction band edge, m* is the effective mass of the electrons, and (3 — {ksT)'^. The 
current (1) contains only the transmitted electrons, i.e., those with sufficient energy to 
surmount the potential barrier 



E:^{x',x")^E:^ix',x")-E',, 



(2) 



where E^[x', x") is the maximum value of the potential profile Ec{x) in the interval [x', x"], 
and E^ is its overall minimum across the sample. The profile Ec{x) comprises the (equi- 
librium) conduction band edge potential and the external electrostatic potential. 

Prom the ballistic electron current J{x',x"), we construct the thermoballistic current 
J{x) at position x inside the sample by summing up the contributions from all intervals 
for which Xi < x' < x < x" < X2, each weighted with the probability of occurrence 
of the interval. For simplicity, we here take this probability in its one-dimensional form 
exp(— |x" — x'\/l), where I is the mean free path for momentum relaxation (momentum 
relaxation length) . We then have 



J{x) — VeNcC < w{xi, X2\ I) 



/•^2 dx' , , 

/ —w[xi,x \ l) 

Jx I 



./3m(^') 



+ 



+ f ^w{x',X2;l) 

Jxi I 

dx' /-^a dx" 



dx f 

Jxi I Jx 



I 



w{x', x"; I) 



(3) 



where 



w{x\ x"\ I) — e~ 



-x'\/l^-l3E^ix',x") 



(4) 



The quantities //i_2 = A*(2^i,2) are the chemical potentials on the contact side of the in- 
terfaces, i.e., immediately outside of the sample. The thermoballistic current J{x) is not 
conserved by itself; however, when averaged over the sample, it yields [14, 15] the physical 
electron current J, 



X2 — Xi Jxi 

Furthermore, the thermoballistic current entering at one end of the sample must equal the 
current leaving at the other end [14] , 



Substituting expression (3) in condition (5), we can derive Volterra-type integral equations 
for two auxiliary functions Xii^) and X2{x), which are distinguished by discontinuities at 
X = Xi and x = X2, respectively. Using condition (6), wc determine from Xi{^) and X2{x) 
a unique chemical potential ij,{x). It exhibits discontinuities at both interfaces, which are 
related to the Sharvin interface resistance [16]. Prom //(x), the thermoballistic current 
J{x) is calculated via Eq. (3). The current- voltage characteristic is obtained for arbitrary 
values of I and arbitrary Ec{x) in terms of a "reduced resistance" x composed of Xi{^'2) 
and X2{xi). The thermoballistic density n{x) of electrons making up the current J{x) is 
constructed in a similar way. 

For electron transport in a homogeneous semiconductor driven by a (constant) electric 
field of strength S, the thermoballistic current J{x) and density n{x) can be expressed 
essentially in terms of the dimensionless parameters x/S, l/S, and eS, where e = l3e\S\. In 
Fig. 2, we show the dependence of n{x) on x/S for eS = 1 and various values of l/S. The 
initial decrease of n{x) with increasing x/S is a ballistic effect that reflects the increase 
of the electron velocity in the electric field and becomes more pronounced as l/S becomes 
larger. The corresponding thermoballistic current J{x) turns out to differ only insignifi- 
cantly from the physical current J. 

3. Spin-polarized transport in semiconductors 

We now extend the unified description by allowing spin relaxation to take place during 
the electron motion across the ballistic intervals. The off-equilibrium spin-polarized current 
J_ {x) = J| (x) — J I {x) is connected with the off-equilibrium spin-polarized density n_ (x) — 




(5) 



J{xt) = J(xj) . 



(6) 
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ni{x) — ni{x) through the balance equation 

dJ-{x) ^ ^ ^ 

dx Ts ' 

where is the spin relaxation time. Using this equation in the ballistic transport regime, in 
which J-{x) — 2ven_{x), we construct the off-equilibrium ballistic spin-polarized current 
J-{x',x";x) at position x in the ballistic interval obtaining 

J_(x', x"; x) = VeN, e-"^^^^''^"^ [A{x') e'^^'^'^/^^ - A{x") e-(-"--)/^»] (8) 

{xi < x' < X < x" < X2), with the ballistic spin relaxation length Ig = 2veTs. The 
parameter Is comprises the overall effect of the underlying microscopic spin relaxation 
mechanisms. The function A{x) is the "spin transport function" defined as 

A{x') = e-'^[^°-'^(^')l a^x') ; (9) 

here, a-{x') — a-[{x') —ai{x') is the off-equilibrium "spin fraction excess" at the equilibra- 
tion point x', which is defined in terms of the spin fractions aii{x'), with ai{x')+ai{x') — 1. 
The spin fractions are related to the spin-resolved thermoballistic chemical potentials 
l^ui^') via 

= eM-')c,^^(^') . (10) 

Proceeding as in the spinless case, we now sum the (weighted) contributions of the ballistic 
spin- polarized current (8) over all randomly distributed intervals [x',x"]. The resulting 
expression for the off- equilibrium thermoballistic spin-polarized current J-{x) is of the form 
(3), but with the terms in brackets therein replaced with those obtained by evaluating the 
bracketed term in Eq. (8) for the different cases. A similar expression is found for the 
off-equilibrium thermoballistic spin-polarized density n-{x). 

Substituting the thermoballistic expressions for J_(x) and n^{x) in Eq. (7), we arrive 
at a linear, Fredholm-type integral equation for the spin transport function A{x), 

W{xux;lJs)Ai + W{x,X2;l,ls)A2 

- Wo{x^, X2; x; I) A{x) + / — W{x', x; I, Q A{x') = , (11) 

where 

W{x', x"; I, Q = w{x', x"; I) e"!^''-^'!/'^ , (12) 



l'X2 

Wo(xi,X2;x;l) ^ w(xi,x;l) +w{x,X2;l) + —■w{x',x;l), (13) 

J XI I 

and 2 = A{xi^2)- The solution of the fundamental equation (11) for xi < x < X2 
determines the spin-polarized electron transport inside the semiconducting sample, and is 
obtained in terms of the values of Ai and A2 on the contact side of the interfaces at the 
ends of the sample. The latter are given by the spin fraction excesses ai^2 — (2^1,2) and 
the chemical potentials /xi 2 = /u(a;i 2) in the contacts via Eq. (9). 

For field-driven transport in a homogeneous semiconductor, Eq. (11) can be converted, 
by twofold differentiation and elimination of the quantities Ai and A2, into a homogeneous 
integrodifferential equation. Within a judicious approximation, the latter equation can be 
reduced to a generalized spin drift-diffusion equation of the form 

^ + ^^^^^ ^ + ^^^^^ ^^^^ = ° ' ^^^^ 
with coefficient functions bi{x) depending linearly on the function exp[— (e + l/l){x — xi)]. 
In the diffusive regime l/lg <^ 1, where bo{x) = 1, bi{x) = e, 62(2^) = —^/{^h), Eq- (14) 
reduces to the standard spin drift-diffusion equation [12] if there the intrinsic spin diffusion 
length L is identified with 

With A{x) determined by solving the integral equation (11) or, in special cases, sim- 
plified equations like Eq. (14), we can evaluate the thcrmoballistic spin-polarized current 
J_(x) from the analogue of expression (3). Dividing by the (total) thcrmoballistic current 
J{x) given by Eq. (3), we obtain the current spin polarization 

W = ^ . (15) 
The density spin polarization is calculated analogously. 

4. Ferromagnet-semiconductor heterostructures 

We now consider spin-polarized transport in heterostructures formed of a homogeneous, 
nonmagnetic semiconductor and two ferromagnetic contacts, which are treated as fully 
degenerate Fermi systems. Equating the splitting ji_{x) = /i|(x) — /U|(x) of the spin-up 
and spin-down chemical potentials in the ferromagnets at a; = 0:1^2, respectively, with that 
of the semiconductor on the contact side of the interfaces, we find, using Eq. (10), 

[/^-(2^1,2)]ferromagnet = [A*- (2^1,2)]semiconductor = i ~^ ~] • (16) 

P \i-Q;i,2/ 
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This condition allows the current spin polarizations Pj{xi^2) on the ferromagnet sides of 
the interfaces to be expressed through the spin fraction excesses ai^2- Neglecting spin- 
flip scattering at the interfaces, we invoke continuity of the polarization at the interfaces. 
We then obtain two coupled nonlinear equations from which ai and a2, and hence J_(x) 
and Pj{x), can be calculated in terms of the material parameters characterizing the fer- 
romagnets and the semiconductor, the bulk polarizations Pi and P2 in the left and right 
ferromagnet, respectively, and the physical current J. Spin-selective interface resistances 
can be introduced via discontinuities of the chemical-potential splitting (x) at the inter- 
faces [8, 10, 12]. 

In the case of zero bias, e — 0, the solutions of Eq. (14) are A{x) oc exp(±a;/L), and 
the position dependence of the current spin polarization inside the semiconductor is given 

by 

P_,(^) ^ ?!^' \c, e-(---i)/^ - C2 e-(^-^)/^l , (17) 
ijj'- 

where L — ^JWs is the generalized spin diffusion length (or "polarization decay length" ) , and 
I — lis I {I + h)- The coefficients Ci^2 are simple functions of the spin fraction excesses ai^2 
which, in turn, are given explicitly in terms of the material parameters, the polarizations 
Pi,2, and the current J. 

In Fig. 3, the zero-bias current spin polarization Pj{x) for a symmetric ferromagnet- 
semiconductor-ferromagnet heterostructure with S — 1 /im at T — 300 K is shown as a 
function of x for various values of the momentum relaxation length I and for zero as well as 
nonzero interface resistances. For the parameters of the ferromagnets, the values 10^ 
cm~^ for the bulk conductivities, 60 nm for the spin diffusion lengths, and 0.5 for the bulk 
polarizations Pi,2 have been adopted from Ref. [12]. With a look at recent experiments on 
the spin dynamics in n-doped GaAs [17], the values 0.067 rUe for the effective electron mass 
m*, 2.0 X 10^^ cm~^ for the equilibrium electron density, and 1 //m for the ballistic spin 
relaxation length Ig have been chosen for the material parameters of the semiconductor. 
We are aware of the fact that, by considering a specific semiconducting material with fixed 
doping concentration, one essentially fixes the value of the momentum relaxation length /. 
Therefore, when varying I in a fairly broad range, we assume the above parameter values 
(or, at least, their order of magnitude) to be representative for a class of semiconducting 
materials that differ in the strength of the impurity scattering and hence in the magnitude 
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of L 

The momentum relaxation length I affects the results shown in Fig. 3 in a twofold 
way. (i) It determines the conduction in the semiconductor and thus the conductivity 
mismatch with the ferromagnets. For small values of I, this mismatch is large, leading to 
a small injected current spin polarization Pj{0). (ii) It determines the polarization decay 
length L, so that for small I the polarization dies out rapidly inside the semiconductor. A 
substantial degree of polarization all along the semiconductor is achieved when the value 
of I is increased up to a length of the order of the sample length, in which case the ballistic 
component becomes prevalent. Figure 3 also shows that, by introducing appropriately 
chosen spin-selective interface resistances, one may offset the suppression of the injected 
polarization due to the conductivity mismatch for small I; however, the rapid decay of the 
polarization inside the semiconductor cannot be prevented in this way. 

For nonzero, constant electric field and -S/L — > oo, i.e., disregarding the effect of the 
right ferromagnet, we have for the current spin polarization Pj{xi) injected at the left 
interface 

Pj{xi)^Xrjai, (18) 

where x is the reduced resistance entering the current- voltage characteristic, and the quan- 
tity Tj involves the function A{x), which is obtained by numerically solving Eq. (14). Since 
only the interface at x — Xi enters into consideration, a single nonlinear equation has to 
be solved to determine the spin fraction excess ai as a function of the polarization Pi. 

In Fig. 4, we show the dependence of Pj{xi) on the electric-field parameter e for vari- 
ous values of the momentum relaxation length I, the remaining parameter values being the 
same as in Fig. 3. In conformity with the drift-diffusion results of Ref. [12], the injected 
polarization generally rises with increasing e; however, as in Fig. 3, the main effect is due 
to the variation of I. 

5. Concluding remarks 

We have outlined the principal ideas of a theory of spin-polarized electron transport in 
ferromagnet-semiconductor heterostructures. It generalizes previous theoretical treatments 
based on the drift-diffusion model by introducing the momentum relaxation length in the 
semiconductor as a new degree of freedom, thus allowing a systematic study of the interplay 
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of spin relaxation and transport mechanism. By considering illustrative examples, we 
have shown that the momentum relaxation length has a significant infiuence both on the 
polarization injected at a ferromagnet-semiconductor interface and on the decay of the 
polarization inside the semiconductor. To study in detail the influence of the transport 
mechanism on spin-polarized transport (in particular, when ballistic effects take over), the 
identification and design of classes of novel semiconducting materials is called for. In this 
way, new possibilities to improve the efficiency of spintronic devices may open up. 
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Figure Captions 



Figure 1: Schematic diagram showing a semiconducting sample of length S enclosed be- 
tween two plane-parallel ferromagnetic contacts. Illustrated are expression (3) for the 
thermoballistic current J{x) and its analogue for the thermoballistic spin-polarized 
current J_(a;). 

Figure 2: The thermoballistic density n{x) inside a homogeneous semiconductor for con- 
stant electric field, plotted versus x/S (assuming xi — 0) for eS = 1 and the indicated 

values of l/S. The density is normalized to the constant value it assumes in the dif- 
fusive limit l/S 0. 

Figure 3: The zero bias (e 0) current spin polarization Pj{x) (assuming xi = 0) along a 
symmetric ferromagnet-semiconductor-ferromagnet structure with S = 1 /xm for the 
indicated values of the momentum relaxation length I. The solid curves correspond 
to zero interface resistance, the dashed curves to interface resistances of 10~^ Q cm^ 
for spin-up electrons and 2 x 10~^ Q cm^ for spin-down electrons, respectively. For 
the remaining parameter values, see text. 

Figure 4: The injected current spin polarization Pj{xi) for 5'/L — > oo as a function of 
the electric-field parameter e for the indicated values of the momentum relaxation 
length /. The solid curves correspond to zero interface resistance, the dashed curves 
to interface resistances of 10~^ fl cm^ for spin-up electrons and 2 x 10~^ fl cm^ for 
spin-down electrons, respectively. For the remaining parameter values, see text. 
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